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ABSTRACT 

For radio pulsars in orbit with a compact companion, pulsar timing observations have proved 
to be a powerful tool for identifying the physical nature of the companion. Unfortunately, perhaps 
the most intriguing system where such a tool could be used, a pulsar in orbit with a black hole, 
has yet to be discovered. 

In this paper we give a detailed investigation of what one can learn about the black hole 
companion via timing observations of the pulsar. We present an analytic calculation for the 
propagation delay caused by the frame-dragging effect and show that it has the same functional 
behavior as the modulation of the observed rotational phase of the pulsar caused by the deflection 
of the radio signals in the gravitational field of the companion (bending delay). Thus, contrary 
to statements of other authors, the frame-dragging delay is unlikely to be separately measurable 
in pulsar binaries where the companion is a stellar mass black hole. We demonstrate, however, 
that the precession of the binary orbit caused by the relativistic spin-orbit coupling can lead to 
observable effects which can be used to set a lower limit to the black-hole spin, or possibly allow 
the determination of its magnitude and orientation. We give parameter estimates for two possible 
systems, a lOM© black hole in orbit either with a young (~ 0.1 s) pulsar or with a millisecond 
pulsar. Finally, we discuss the measurability of the quadrupolc moment of the rotating black hole 
companion which would test the presence of a Kerr black hole. As an interesting side result of 
our calculations, we can give a further argument why the companion of PSR J0045-7319 cannot 
be a Kerr black hole. 



Subject headings: Black hole physics — gravitation — binaries: general — pulsars: general 
pulsars: individual (J0045-7319) 



1. Introduction 

So far the best arguments for the existence of stellar mass black holes (BHs) are based on dynamical 
mass estimates in X-ray binaries. The measurement of absorption-line velocities of the secondary star allows 
us to determine a lower limit to the mass of the compact companion. If the mass of the companion exceeds 
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the calculated maximum mass of a neutron star (~ 3M Q ) we call it a BH candidate (see Wijers 1996 for a list 
of BH candidates). But to argue convincingly that the companion is indeed a BH one must rule out possible 
alternatives to stellar mass BHs. In case of cold, non-rotating neutron stars (NSs) the permissible maximum 
mass depends on a rather complicated equation of state (EOS) for condensed matter at and above nuclear 
density. Since our knowledge of the EOS for matter above nuclear density is yet imperfect the maximum 
masses of NSs range widely, from 1.5M Q to about 2.5M Q (Friedman, Ipscr, & Parker; Komatsu, Erguchi 



fc Hachisu 1989| ; Lattimer et al. 199C; Weber & Glendcnning 1992). This is well below the dynamical mass 
estimations for many of the suspected BH binaries. However, if we are very conservative and assuming 
complete ignorance about the EOS above a density of 10 14 g cm~ 3 , dropping the causality condition (i.e. 
allowing a dispersive medium), and allowing rotation, then the maximum mass goes up to about 14M© 



( |Fricdman fe Ipscr 1987| ). A compact object with a maximum mass of 14Af Q can explain all the dynamical 
mass estimations in BH candidates. There are other ways of exceeding the dynamical mass limits, such as 



rather exotic compact stars like Q-stars (Bahcall, Lynn, & Selipsky; Miller, Shahbaz, & Nolan 1998) or boson 
stars made from a strongly self-interacting scalar field ( |Colpi, Shapiro, fc Wassermaii 1986| ), the latter having 
in principle no maximum-mass limit. Finally, one might even think of abandoning general relativity as the 
correct theory for the description of the strong gravitational fields inside a neutron star (although pulsar 
timing tests of the strong equivalence principle support the validity of general relativity in these strong field 



regime (Damour fc Schafer 1991; Wcx 1997) 



Recent advances in the theory of accretion physics (e.g., the advection-dominated accretion flow model) 
lead to the conclusion that the matter accreting onto the compact companion in some of the BH candidates 
does not hit a hard surface. This may indirectly indicate the presence of an event horizon through which 



matter and radiation can fall in but from which nothing escapes, a fundamental property of BHs (Narayan 



Garcia, & McClintock 1997). At present the concept of an event horizon in these BH candidates has rather 
the character of a "most natural explanation" than of an inevitable necessity. Certainly the absence of a 
hard surface does not prove the presence of a BH as predicted by Einstein's theory of gravity, i.e. a Kerr BH. 

Astrophysical BHs are expected to rotate. The spin of the BH gives rise to a so called gravitomagnctic 



field (Thorne & Hartle 1985) which causes the relativistic dragging of inertial frames in the vicinity of the 
BH. The dragging of inertial frames has two important consequences for an accretion disc around the BH. 
First, the radius of the last (marginally) stable orbit of the accretion disk is a function of the BH spin 



(Bardccn, Press, & Teukolsky 1972). Secondly, any deviation of particle motion from the equatorial plane 
will cause a precession of the particle orbit ( Lense fc Thirring 1918 ; Bardeen fc Pettcrson 1975 ; Shakura 



1987). In principle, these effects allow the determination of the spin of the accreting BH using observations 



of low-energy X-ray spectra (Zhang, Cui & Chen 1997) or high frequency quasi-periodic oscillations (Cui 



Zhang, & Chen 1998). While it is still unclear whether the observed quasi-periodic oscillations are a result of 
the Lense- Thirring precession of inner parts of a warped accretion disk, spectral analysis seems to show the 
presence of the gravitomagnetic field of the BH ( Markovic fc Lamb 1998 ; McClintock 1998j ). Nevertheless, 



the arguments are still require certain plausible assumptions (e.g., the accretion disk reaches all the way 
down to the last stable orbit, etc.) and the precision with which the spin is determined is rather poor, even 
in systems where all the input parameters are comparably well known. 



The discovery ( Hulse fc Taylor 1975|) and continuous observation of binary pulsars (i.e. pulsars in orbit 
with a gravitationally bound companion) has opened an entirely new testing ground for gravity theories. So 
far binary pulsars provide the only laboratories where one can probe the gravitational radiation properties 
and strong field aspects of relativistic gravity ( Taylor fc Weisbcrg 1989) ; Damour & Taylor 1992; Damour & 
Esposito-Farese 1996a; Damour & Esposito-Farese 1996b). A useful review on pulsar timing and relativistic 



- 3 - 



gravity can be found in |Taylor 1993 , where the reader can find the basic ideas of this subject. 



A pulsar orbiting a BH in a close orbit would certainly be a high precision laboratory for BH physics 



(Paczynski & Trimble 1979). After a short observation period it would be possible to measure the relativistic 
advance of periastron, oj, and thus obtain the total mass of the binary system to high precision. At this stage 
the mass of the BH would be fairly well constrained since the mass of the pulsar is unlikely to exceed 2Mq. 
The measurement of a second post-Keplerian (PK) parameter would then yield the mass of the pulsar, m p 



(= A4 p Mq), and the mass of the companion, m c (= A4 c Mq). It was even suggested ( Narayan, Piran, fc 



Shcmi 1991 , Laguna fc Wolszczan 1997 ) that a pulsar in a close, nearly edge-on orbit, with a spinning BH 
might be used to probe the rotation of the BH via the timing delay produced by the dragging of inertial 
frames caused by the spin of the BH (frame-dragging propagation effect). Unfortunately so far none of the 
binary pulsars seems to orbit a stellar mass BH of typically J> 10M Q (see Fig. |l|). 

However, it is possible that pulsars in close orbit with BHs have already been missed in previous pulsar 
surveys due to the Doppler smearing of the signal during the integration. Whilst it is true that this effect 



can be largely removed by the application of an "acceleration code" (Johnston & Kulkarni 1991), only a 



number of globular-cluster searches have utilized this technique due to previous limitations in computing 



power (Anderson 1992). The continual improvement and availability of more powerful computers, however, 



means that many more surveys are now able to incorporate acceleration searches (Lorimer 1998) 



In this paper we discuss what one can learn about a stellar mass BH by timing its pulsar companion. 
We focus, in particular, on the observational determination of the BH's rotation (angular velocity and 
orientation) and assume throughout the paper that in addition to the five Keplerian parameters, at least 
two PK parameters are measured with high precision, so that we know A4 P , A4 C , and therefore sini, i.e. 
the angle of orbital inclination, i, up to the ambiguity i — > n — i. We begin with a brief introduction to 
the definition of the spin and quadrupole moments of a rotating body in general relativity and relations 
between them (§2). In §3 we give an analytic treatment of the frame-dragging propagation effect and 
show that it is practically impossible to observe this effect in stellar mass BH-pulsar binaries. In §4 we 
concentrate on the orbital dynamics of a pulsar orbiting a rotating BH and give a detailed investigation 
of the precession of the pulsar orbit caused by the spin (relativistic spin-orbit coupling) and quadrupole 
moment (classical spin-orbit coupling) of the BH. In §5 we calculate the secular changes in observable 
quantities caused by the relativistic spin-orbit coupling. We show that, in principle, one can extract the spin 
of the BH from nonlincar-in-time changes of the observables which can be approximated by polynomials in 
time with sufficient accuracy. We discuss additional effects that can cause secular changes in the observed 
parameters reducing the measurement accuracy for the BH spin. In §6 we outline the prospects of BH spin 
determination once a BH-pulsar binary is discovered. In §7 we present a method of extracting the quadrupole 
moment of the rotating BH from the orbital dynamics of the pulsar and evaluate its actual measurability. 
We summarize in 58. 



2. Spin and quadrupole moment of a compact body 



The external metric of a stationary, axially symmetric body can be written in the standard form (Bardeen 
fc Wagoner 197l[) 



ds 2 = _ e 2- c 2 dt 2 



R 2 sin 2 eB 2 e^(d(j) - uj dt) 2 + e 2a {dR 2 + R 2 d9 2 ) , 



(1) 



where c is the speed of light; R, <f>, 6 are ordinary spherical coordinates, and the potentials v, B, uj, and a 
are functions of R and 6. Butterworth and Ipser (1976) have calculated the asymptotic behavior of these 
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potentials for large R. In particular, the asymptotic behavior of the potentials v and u) is (using c = G = 1) 

(2) 
(3) 



M 1 fM\ 3 fM\ 3 , ns „ ( 1 



i? 



R 4 



(Butterworth & Ipser 1976, Laarakkers & Poisson 1997); M is the body's mass and \ is a dimensionless 
measure of the body's angular momentum (spin), S: 



X 



c_ _S_ 

G AP 



(4) 



The dimensionless parameter q in equation (Q) is related to the quadrupole moment, Q, of the rotating body 

c 4 Q 



<1 



G 2 M 3 



In the Newtonian limit one has 



Q = J Q(R',6')R' 2 P 2 (cos9')dV , 



(5) 
(6) 



where £> is the mass density of matter inside the star. The Newtonian theory of self-gravitating, rotating 
bodies predicts a certain relationship between the angular velocity of the body's rotation and its quadrupole 
moment due to the oblateness in the mass distribution ( |Chandrasckhar 1969 ). For the same reason, one 
now may expect a complicated relation between the rotational parameter of a star, \t and its quadrupole 
moment, g, in relativistic gravity theories. For a rotating Kerr BH in general relativity X < 1 for if one 



had x > 1 it would give an unacceptable naked singularity (see Hawking fc Ellis 1973 ). Furthermore, for a 
Kerr BH the quadrupole parameter, g, is uniquely determined by the rotational parameter x ( Thome 198C| , 
Thornc, Price fc Macdonald 1986Q : 

g = -x 2 - (7) 

Similarly for solid extended bodies, Laarakkers and Poisson (1997) have computed the quadrupole 
moment of rotating neutron stars in the interval between 1.0 and 1.8 solar masses using four different EOS. 
They have found the maximum values for x lying in between 0.62 (softest EOS) and 0.73 (stiffest EOS). 
For the quadrupole moment of the mass distribution, g, they have derived a dependence on x which is well 
reproduced by the relation 

q ~ -Cx 2 , (8) 
where the constant C takes a value between 2.0 and 12.1 depending on the mass of the NS and the EOS. 
For spinning boson stars with large self-interaction the maximum value for the parameter x seems to 



be between 3 and 4, and the relation between x and q is rather complicated once x exceeds 0.2 (Ryan 1997) 
The following inequality, however, always holds: 



q < -i0x z 



(9) 



Comparing Eqs. (0) - (^|) one sees that simultaneous measurements of the mass, M, spin parameter, Xj 
and the quadrupole moment, g, of the pulsar's compact companion can lead to a unequivocal identification 
of its physical nature]^] Such measurements seem to be a straightforward way to observationally verify the 
existence of a rotating Kerr BH. 



5 We are not aware of any calculations done for Q-stars. 
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3. The frame-dragging propagation effect in highly inclined binary systems 

In Newtonian gravity the rotation of a star contributes only indirectly to its gravitational field through 
the rotationally induced oblateness of the mass distribution. The centrifugal flattening of the surfaces of 
equal mass density caused by the rotation of the star gives rise to a gravitational quadrupole field. In general 
relativity, however, not only mass itself contributes to the gravitational field. The rotation of a star gives 
rise to a mass current which contributes in a specific way to its gravitational field. The rotation is the source 
of the gravitomagnetic field of the rotating body ( Thome, Price fc Macdonald 1986| ). This gravitomagnetic 



field influences the motion of test particles and the propagation of light in the vicinity of the rotating body 



see, for instance Ciufolini & Wheeler 1995, Mashhoon 1997 and references therein). It is well-known that 



BHs, although purely geometrical in their nature (i.e. a solution of Einstein's vacuum-field equations), are 
physical objects and as such also can carry angular momentum which gives rise to a gravitomagnetic field. 
In fact, the strength and configuration of the gravitomagnetic field of the BH is used to define its rotation. 



The rotation of the BH will influence the propagation of photons in a specific way ( Bardeen 1973 ). For 
pulsar-timing experiments its influence on the propagation time of radio signals are of particular interest. 
We call this contribution to the propagation time of the pulsar signals the frame- dragging propagation effect. 

Recently, Laguna and Wolszczan (1997) have suggested on the basis of numerical ray-tracing calculations 
that pulsar-timing experiments could measure the gravitomagnetic field of a rotating BH companion by its 
influence on the propagation of the pulsar signals emitted from different points of the pulsar's orbit while 
the pulsar is passing through the point of the superior conjunction with the companion. In this section we 
give a more detailed investigation of this frame-dragging propagation effect based on analytic calculations 
and show the problems preventing its direct measurement. 

For a discussion of the effect in question it is convenient to use the metric of Kerr spacetime, written 



in terms of Boyer-Lindquist time t and asymptotically Cartesian spatial coordinates X J . One has (Thornc 



Price fc Macdonald 1986| ): 

ds 2 = -a 2 c 2 dt 2 + g jk {dX j + f3 j dt)(dX k + (3 k dt) , (10) 

where at large distances from the BH one finds for the lapse, a, and shift, (3 l , functions the following 
expressions (Thorne, Price fc Macdonald 1986| ): 



where R = |X| is the distance from the BH to the current point in space. The shift function (3 represents 
the vector potential of the gravitomagnetic vector field which drags local inertial frames into rotation (frame 
dragging effect) . Since for photons ds 2 = the contribution of the frame dragging effect to the propagation 
time of a photon from the pulsar to the observer is 

[3-dX 2G (S x X) ■ dX 
dt FD - — g- = - 3 . (12) 



After substituting into the right-hand side of equation (|12| ) the unperturbed trajectory of the photon, which 
is a straight line, the equation is easily integrated in a coordinate system where X = (x, 0, b). The constant b 
is the minimum distance of the unperturbed trajectory of the photon to the BH^j Consequently, the timing 



Details abou t the propagation of light rays in the stationary field of relativistic gravitational multipolcs can be found in 



(Kopeikin 1997e) 
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delay caused by the frame-dragging effect assumes the form 



2G a l f°° dx 2G S V (■> X 



AFD . __ s „^ = __ f ^ _ , (13 ) 

where ro = |Xo| is the distance between the BH and the point of emission of the photon which is located 
at the pulsar's orbital position. We have assumed that the observer is far away from the rotating BH; for 
this reason, the upper limit of the integral in ( |l3| ) was taken to be infinity. Indeed, from equation ([l2]) it 
is already clear that most of the contribution to this frame- dragging (FD) propagation effect comes from a 
region around the BH with a few impact parameters, b, in radius. 

Let R denote the vector pointing from the BH to the pulsar and Ko denote the unit vector pointing 
from the observer to the pulsar (see Fig. 0). Then, R = \tq\, and one has (cf. eq. (43) for I = 1 in Kopcikin 



1997a| ) 

_ 2G S-(KqxR) 
AfD ~ ~? R(R Ko • R) ■ (U) 
Using the angles as defined in Fig. || we find for the FD propagation effect in a binary pulsar system 

Apr, = A~ 1 ^Afd cos i sin [ui + A e (u)] + Bfd cos [u> + A e (u)] X , (15) 

where the function A is defined by 

A = 1 — ecosu — s[sinw(cosu — e) + (1 — e 2 ) 1 / 2 coswsinu] , (16) 
the eccentric anomaly angle is 



A e (u) = 2 arctan 



1 + e x 



1 



tan 



(17) 



and 



2vr\ 2/3 Ml 



A FD = +4T^j { M P + M.)^ XS[nX ' C ° Sfl - (18) 
B FD = -4^(|) 2/3 ( ^^ )1/3 xsinA.sin^. (19) 
When Ai p <C A4, we find the following numerical estimate for the constant factor 

T =« (HW /3 *g „ ,0.000! „, (-*-)"" (»•)'" . ,20) 

To illustrate the expected strength of the FD propagation effect we assume a pulsar in a circular orbit 
where n m = and the inclination angle i is close to 90°, i.e. e = 0, u> — 0, A e — u, s — sini ~ 1, and 
cosi <C 1. Then Afd 7^ 0, £>fd = 0, and 

sin(o; + A e ) 



A 1 — s sin u 



(21) 



This function has maximum at u = n/2 and minimum at u = 3n/2. Thus, for the maximal magnitude of 
the FD effect we have: 

/A , 0.0008ms/ Pb Y 2/3 {M.\ 5/3 . , 
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Hence, for BH binary pulsars with an orbital inclination i very close to 90°, the FD propagation effect 
will have a measurable influence on the TOAs, provided that the pulsar is a millisecond pulsar where one 
can expect (with present-day technology) a timing accuracy of better than one microsecond. This fact was 
already pointed out by Laguna & Wolszczan (1997) using numerical methods to study the propagation of 
the pulsar signals in the spacetime of the Kerr-BH companion. In particular, for a binary pulsar with an 
orbital period Pb about 1 day and an extreme 10M Q Kerr BH companion with \ — 1 the strength of the FD 
propagation effect will be one microsecond if cosi < 0.001 or \i — 90° | < 0.05°. 

A measurement of at least one of the parameters Afd or Bfd will give values for either ^sinA, cos 77. 
or xsin A. sin 7/.. Since we have a priori no idea about the orientation of the BH spin it is not possible to 
extract \ from these measurements. But we get a lower limit on which has to be less than one, otherwise 
it is not a Kerr BH. In the lucky case of having good values for both parameters, Afd and Bfd, we can 
calculate 77. and consequently \ sm ^» which gives a definite lower limit for the \ of the companion. An 
upper limit to \ can be set only with a certain probability excluding values of A, being close to either 0° or 
180°. 

Unfortunately, the measurement of the FD propagation effect is complicated (and may be impossible) 



by a competing effect which occurs during the superior conjunction, the bending delay (Doroshenko & 



Kopeikin 1995|) The delay originates because of the modulation of the pulsar's rotational phase by the effect 



of gravitational deflection of light in the field of the pulsar's companion. The expression for the bending 



delay in the framework of general relativity is given by the formula (Doroshenko & Kopeikin 1995) 

A B = A _1 [^l s cosisin(w + A e ) + B b cos(lu + A e )] , (23) 

where 

T| /3 /27r\ 2/3 M. cosTTp 



AB + itv p \P b ) (Mp + M.) 1 / 3 sinA p ' (24) 

= T| /3 [2n\ 2/3 M. sinr ?p 
3 nu p \P b J (Mp + M.) 1 / 3 sinA p " 1 ' 

Here, v p is the rotational frequency of the pulsar, and the angles X p , r] p have the same meaning as A., 77. but 
define now the orientation of the pulsar's angular velocity vector in space. For M p <C M. we find 

zgvyf ro , ( o.oo 07)ts) f m ( ft r'VAM 2 ' 3 , (26) 

tti/j, \PbJ (Mp + M.) 1 / 3 V P ; \lmsj Vldayy \ 10 / y ' 

where P denotes the rotational period of the pulsar. 

The bending delay has the same functional dependence on the orbital motion as the FD propagation 
effect and, for this reason, one can observe only a linear combination of the FD effect and the bending delay. 
Hence, in timing observations the parameters of the FD effect are not separated from those of the bending 
delay. For pulsars with rotational periods close to 1 ms the amplitudes of these two effects are comparable if 
M. ~ 20. For smaller BH masses and/or longer pulsar periods the bending delay is clearly the dominating 
effect. Only for BH companions with more than ~ 100M Q the FD propagation effect will dominate over the 
bending delay. 



In principle, one may extract r] p and A p from additional pulse structure analysis (Damour & Taylor 



1992). In this case the bending delay can be predicted and subtracted from the total effect. In practice, 



pulse structure observations and subsequent analysis are connected with rather large uncertainties in the 
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corresponding angles. Therefore, while using results of such an analysis one can still hope to get a comparably 
good measurement for A p , the angle rj p affects the polarization pattern of the pulsar only at the negligible 
small x/Pb level (Damour fc Taylor 1992) and, thus, will remain unobservable. 



4. Spin-orbit coupling and precession of the binary orbit 



Eighty years ago Lense & Thirring (1918) pointed out that the gravitomagnetic field of a central rotating 
body will cause a precession of the orbit of a test particle. In the same way, the rotation of one or both 



components of a binary system will cause a precession of the binary orbit (Brumberg 1991). In this and the 
following sections we will show that the observation of such a precession can lead to the direct determination 
of the spin of a BH companion. 



The typical moment of inertia for a pulsar is of order 10 45 cm 2 g ( Arnett fc Browns 1977 ). Thus, 
assuming the rotational period of the pulsar to be just 1 ms we obtain a spin of order 6 x 10 48 cm 2 g/s. Using 
equation we find that the spin of a lOAf© extreme Kerr BH (\ — 1) has a value more than 100 times 
bigger. Since for even a very soft EOS the rotational period of a pulsar should lie above ~0.5 ms (Weber & 



Glcndcnning 1992), the spin of the BH will dominate the orbital precession. For this reason, we neglect the 
spin of the pulsar in our subsequent calculations. The secular precession of the orbit is then given by two 



vector differential equations of first order (Brumberg 1991) 



^prcc X L 



^prcc X A 



(27) 



where L is the orbital angular momentum and A is the Laplace-Runge-Lenz vector which points to the 
instantaneous position of periastron of the precessing orbit. The angular velocity vector f2 proc is the sum of 
the well-known relativistic (post-Newtonian) periastron advance, fipjv, the gravitomagnetic Lense-Thirring 
precession, fig, caused by the coupling of the orbital angular momentum vector to the spin of the BH, and 
the classical precession, Oq, due to the Newtonian coupling of the orbital angular momentum vector to the 
quadrupole moment of the rotating BH 



fiprec = flpN + fls + 



where 



ttpN 



QpN > 

n* s [3(1 



S)L 



% [{5(L ■ S) 2 - 1}L - 2(L • S)S] , 



(28) 

(29) 
(30) 
(31) 



and a hat on a vector indicates the unit vector in the same direction as that of the vector itself. It is 
important to stress that fig and Hq are not equal to the absolute values of the corresponding vectors. (We 
have deliberately used the asterisk to avoid such a possible confusion.) One finds (see Barker & O'Conncll 



1975) 



^2/3 



2tt\ 5/3 



M.) 



2/3 



X T Q 



(1 - e 2 ) 3 / 2 
(1 - e 2 ) 2 



27ry M.(3M P + 4M.) 

2^ 7/3 
Pb. 



2(M P + M.) 

3Ml 
4(M P + M.) 2 / 3 



(32) 
(33) 
(34) 
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It is worth emphasizing that both vectors L and S precess around the conserved total angular momentum 
of the binary system J = L + S, while their absolute values, L = |L| and S = |S|, are conserved quantities 
(averaged over one orbital revolution). 

The precession of the orbital plane and the longitude of periastron are best described by the angles $ 
and \I> defined with respect to the invariable plane, i.e. the plane perpendicular to J (see Fig. ||for details). 
From the identities 



Jlp rcc xL = L = $JxL. 

n prec x A = A = ($J + \j/L) x A 

one obtains the precession in terms of the angles $ and 

$ = &s + $Q = const. , 



where 



-n* s 



sin 9 
sin 9 j 
sin 29 
sin 9 j 



and 



where 



* = *pjv + *S + *Q = const. 



^PN — — ^GRO > 

*s = ilg (2 cos + sin cot 0j) , 

* Q = fi^ (| + § cos 20 + sin 20 cot 0j 



(35) 
(36) 



(37) 

(38) 
(39) 

(40) 

(41) 
(42) 
(43) 



In particular, eqs. 
19951) . 



(|43|) can be found in Smarr & Blandford (1976) (see also Lai, Bildsten, fc Kaspi 



It is worth noting that the rate of the precession depends on two angles 9 and 9j. However, given the 
spin of the BH, S, one can express 9j as a function of 9, since Lsin9j = Ssin(0 — 0j). We find 



1 



1 



sin 9j 



L 2 L 
-^+2 - cos0 



1/2 



(44) 



cot0j = cot 9+——, (45) 
b sin0 



where 

5_ x4 /3 /2^\ 1/3 A^.(M P + M.) 1/3 



L (i_ e 2)i/2 ^ Pb y Mp 

Hence, the rate of the orbital precession may be expressed as a function of only one angle. 



(46) 



We conclude this section by giving numerical values for the various precessional effects. Assuming a 
pulsar mass of 1.4M and a Kerr BH mass A4 t J> 10 allows us to approximate equations (]32|) to ([54]) by the 



-10- 



expressions 



,-3 X ( Pb \ 2 (M, 



°S " °' 9 x 1(r J (rf^ [ufc) Uf J des/yr < (48) 

- 0.5 x 10-- - (rS^ ( ^ deg/yr . (49) 



Now it is easy to see that the precession in ^ is clearly dominated by Q PN while that in $ is dominated by 
the relativistic spin-orbit coupling fig. The precession Hq caused by the quadrupole moment of the BH is 
three orders of magnitude smaller than the relativistic spin-orbit precession and, therefore, will be omitted 
from the following discussion. 



5. Spin-orbit coupling and observable quantities 

The angles $ and 'J', which change linearly in time, are not directly observable in analyzing pulse arrival 
times. Instead, one can extract from the timing observations the projected semi-major axis of the pulsar 
orbit, x = a p sini/c, and the longitude of periastron, lo, which are connected with $ and \& through the 
trigonometric relationships (see Fig. ^ 

cosi = cosij cos#j — sin#j s'mij cos $ , (50) 

and 

sinzsinw = (sin#j cosij + cos 6j sin ij cos $) sin ^ + sinij sin $ cos \& , (51) 

sin i cos lo — (sin 0j cosij + cos 6j sin ij cos $) cos *f> — sinij sin $ sin \& . (52) 

Consequently, due to the non-linear character of these relationships, the lincar-in-time precession of <3? and 
'J will cause a non- linear- in-time precessional evolution of the observed parameters x and lo. Although the 
changes in lo and x can be expressed in a closed analytic form, a Taylor expansion in powers of time t — To, 
where To is the time of periastron passage, is more suitable for observational purposes, as will become clear 
in the next section. Thus, one has: 

hf + ... , (53) 
tof + ■■■ , (54) 
where xq, loq are initial values of x, lo at the initial epoch To, 

±p- c = XS + XQ (55) 
cj PI ' OC = U) GR0 +Lb S +LOQ , (56) 

and subindeces S and Q mean that the corresponding quantity is caused by the influence of the BH spin or 
the BH quadrupole moment, respectively. The quantity logro describes the standard relativistic periastron 
advance. 



x picc (t) = Xo + x P^(t-t ) + -xP ICC (t- 

LO PrCC {t) = LO +LO PrCC {t-t ) + l^ CC (t- 



- 11 - 



As mentioned previously, the contribution of the quadrupole moment of the BH companion to any 
observable parameter is extremely small. For this reason, we again omit, in what follows, all quantities with 
the subindex Q. Further simplification is achieved if one makes use of the fact that, in general, S <C L and, 
thus, as a consequence of equation 9j rs SsinO/L <C 1. Indeed, when M, ;> 10, 

S nm X ( Pb Y 1/3 (M.^ i/3 



Therefore, even for a binary system with an orbital period P& of 0.5 days, an eccentricity e of 0.8, and a 
2OM BH companion 9,j will be smaller than ~ 3°. Taking this into account equations ([30]) and ( |5l| ) assume 
the form: 

sin i — sin i j + 9j cos ij cos $ + O(0j) (58) 

and 

w = + 6j cotijam$ + O(0 2 j) , (59) 



(cf. |Wex 1998| ). Thus, one obtains, to leading order in 9j, 



Co s ^ -Bg x sin 9 cos $ - 2A* S \ cos 6 , (60) 

— ~ -B* s xsin6»sin$ , (61) 

lo s ^ +Cg x sin sin $ , (62) 

^ ~ -C5 xsin6»cos$ , (63) 
where $0 is the numerical value of the angle $ at the initial epoch Tq and the constant coefficients 



A s = n*/ X , B* s =A* s coti, C* s = B* s n* s L/S (64) 

are independent of the angles 9, $0, and the spin parameter x- The values for these coefficients are fixed 
once we know the masses of the binary system. 

At the beginning of the timing project one will just be able to see the linear trend in the precession, 
<j, and some time later x. But after a while, depending on the compactness of the binary system, one will 
start to see 10 and 'x. Since Co is clearly dominated by ojgro the first sign for the existence of spin-orbit 
coupling in the binary system will come from the measurement of x. Under favorable circumstances (see 
equation ( |65| ) and related discussion) this allows us to calculate x sin 9 sin $0 if w e assume that the masses 
of pulsar and BH and, thus, the angle of the orbital inclination, i, are known from the measurement of at 
least two PK parameters. Since |xsin^sin$ < Ix we can S e t a lower limit for the spin of the BH or, 
if we find \x\ > 1, a Kerr BH is ruled out. In fact, it is interesting to note that if one assumes that the 
precession of PSR J0045-7319 (see Kaspi et al. 1996) is caused by a BH companion, the analysis presented 
here would give |x| J> 100 which is certainly not a Kerr BHj^| The measurement of x will then allow the 
separate determination of the angle $0 and x sin 9. If wc can determine the masses of the bodies comprising 
the binary system without measuring Co (e.g., using 7 and Pb) then, in principle, we can extract Cos from 
the total observed value of Co and finally determine x- It is important to see that neither Co nor even higher 
derivatives of x and lo contain a term xcos(9. Thus the only way to fully determine the spin of the BH is 
the separate measurement of Cos- 



7 There is clear observational evidence that the companion of PSR J0045-7319 is a B star (Bell et al. 1995). The precession 
of the orbit is caused by the large quadrupole moment of the B star (Kaspi et al. 1996) 
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Unfortunately, the precession of the orbit due to spin-orbit coupling is not the only effect that gives rise 
to the parameter x in the timing model. Following Damour & Taylor (1992) and Kopeikin (1994, 1996) one 
finds, in fact, four alternative effects causing a secular change of the observed value, x ohs , of the parameter 



dSA 

dt 



D 
D 



(65) 



The first contribution to the right-hand side of equation ( |65| ) is due to the precession of the binary orbit 
discussed at the beginning of this section. The contribution to the right-hand side of equation (|65| ) from the 
second (gw) term is caused by the shrinking of the pulsar orbit due to the emission of gravitational waves 



by the binary system. Given the masses of pulsar and companion this effect can be calculated (Peters 1964) 
and subtracted. The third (pm) contribution is due to the proper motion of the binary system, which causes 
a gradual change in the apparent geometrical orientation of the orbital plane (Kopeikin 1994, 1996; see also 
Arzoumanian et al. 1996) If the pulsar is close to the solar system and/or has a high proper motion this 



contribution can become quite significant, in which case, additional astromctric information will be needed 
for measurement and subtraction of i pm . The forth contribution to the right-hand side of the equation ( |65| ) 
is due to a varying aberration caused by the deSitter precession of the pulsar's spin while the pulsar orbits 
its companion. The pulsar spin precesses around the direction of the orbital angular momentum, L, with 



angular velocity (Barker & O'Connell 1975) 



Qspin 



T 



2/3 



5/3 



M.(4M P + 3M.) 
2(X P + M.) 4/3 



and, consequently, (Damour fc Taylor 1992) 



deA P ^ P pin cot X p sin 2ij p + cot i cos r\ v 
~dT = ~P~ b (1-e 2 ) 1 / 2 sinAp 



(66) 



(67) 



where P is the pulsar's rotational period. A comparison with equation (|6l]) leads to (excluding values of A p 
close to 0°and 180°) 



x s /x 



~p h ~n% 



(1.(106 | — 

1 s 



Pu 



1 day 



-2/3 



M, 

10 



-1/3 



(68) 



Thus, the varying aberration effect will be of the same order as the precession of the orbit only for slow 
pulsars (P ;> 1 s) in very close orbits with a BH companion. But in this case the varying aberration will 



also show up as a change of the eccentricity ( Damour fc Taylor 1992 ) 

obs 



de A 
dt 



(69) 



and can therefore be observed and subtracted. (The change due to gravitational wave damping (gw) can be 
calculated if the masses are known (see Peters 1964[ ).) 



The last term in equation (pq ) is due to a varying Doppler shift, D, caused by a secular change of the 
distance between the solar-system barycenter and the binary-pulsar system. The effect under discussion is 
non-linear in time and can be caused both by the acceleration of gravitational fields in our galaxy and an 



8 Apart from a contribution to x this effect changes also the observed value for the advance of periastron, u) (Kopeikin 1994, 
1996). 
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apparent acceleration due to the proper motion of the binary system in the sky ( [Bhklovskii 197C ; Damour 
fe Taylor 199l|): 



jj — Kg " (^binary system 



^■solar system) H~ 



cd 



(70) 



Here, Vt denotes the velocity of the pulsar being transverse to the line of sight, d is the distance between 
the pulsar and the solar system, and abi nar y system and a so i ar system are galactic accelerations at the locations 
of the solar and binary systems respectively. The galactic acceleration at the solar location is of order 
2 x 10~ 8 cm s -2 which is a typical value for the acceleration fields in the galactic disk ( parlberg fe Innaiien 
1987; Kuijken fc Gilmorc 1989] for a model of the Galactic acceleration field). Hence, unless the pulsar is 



located in the core of a globular cluster where accelerations are much higher, the contribution of galactic 
gravitational fields to x is negligible. The same is true for the Shklovskii term if we assume typical distances 
(;> 1 kpc) and typical pulsar velocities (few 100 km/s). Thus, we conclude that a varying Dopplcr shift should 
not complicate the interpretation of timing observations of BH-pulsar binaries under ordinary circumstances. 

As pointed out above, for a full determination of the rotational parameter \ of the pulsar companion, 
one cannot make use of Co to calculate the masses of pulsar and BH. Thus, one needs a third PK parameter 
for the separate mass determination which will be most likely Pf, (caused by gravitational wave damping). 



Since (Damour fc Taylor 1991) 




obs 




D 
D 



and 



Pb 



gw 



(-1.6 x 10- 18 s- n 



Pb 
1 day 



2/3 



M. 

10 



2/3 



(71) 



(72) 



one has to be cautious about the contribution of the galactic acceleration and Shklovskii's term. Even if we 
understand these effects to about 10% for BH-pulsar binaries with orbital periods greater than a few days, 
a reliable determination of tjjg will be impossible. Exceptions are cases where the system is seen sufficiently 
edge-on so that one can measure the two PK parameters related with the Shapiro delay, r and s, with 
necessary precision. For PSR B1913+16, the limit on the determination of P^" due to galactic accelerations 
has been clearly demonstrated by Damour & Taylor (1991). 



For a BH-pulsar binary the deSitter precession of the pulsar spin, equation ( 66 ) , is approximately lvpn/2. 
For short orbital periods this can be a few degrees per year or even more. Then, if the pulsar spin is tilted 
with respect to the orbital angular momentum the deSitter precession can lead to significant changes in the 



pulse profile (cf. Weisbcrg, Romani, fc Taylor 198S ; Kramer 1998 ) which might, in principle, limit the timing 
accuracy. In extreme cases the pulsar could even cease to beam in our direction after a couple of years. 

Only if the influence of all these additional effects listed above is well understood we can hope to perform 
a spin determination for the BH companion of the pulsar. 



6. Parameter estimations 

So far we have shown that in principle, one can use the pulsar timing observations of a BH-pulsar binary 
system to determine the mass and spin of the rotating BH companion. In the previous section, we discussed 
in detail the extraction of the BH mass and spin from the observed values of Keplerian and PK parameters. 
However, the number of measured parameters one hopes to determine from timing observations depends 
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critically on both the accuracy of the measured TO As and the compactness of the binary system's orbit. For 
example, we show that an ordinary one-second period pulsar in a one-year orbit with a 10 Mq BH will allow 
at most the determination of the mass function and sets only a lower limit on the mass of the BH companion. 
However, for a millisecond pulsar in an eccentric one-day orbit with a 10M Q BH, the masses of the pulsar 
and companion can be determined with high precision after a month of regular timing observations from the 
measurement of two PK parameters (most likely Co and 7). After a one-year observing campaign one can 
expect to see additional relativistic effects, such as gravitational radiation damping and relativistic spin-orbit 
precession. In this section, we explore the potential measurability of various relativistic effects in BH-pulsar 
binaries. We use a standard tool, the information matrix, to asses the errors of our parameter estimation 
method. 

Let N denote the total number of measured TOAs, r, (i = 1,2, ...,N), and M the number of fitting 
parameters, £ a (a = 1, 2, M), in the timing model used for the a priori estimate of A/"(t, £) the number of 
the pulse arriving at the time r. The parameters are determined using a least-square minimization method 
where the goodness of fit parameter, \ 2 1 is given by 



x 2 = E 



N 

' n,. 



(73) 



Here v v is the pulsar's rotational frequency, rij is the closest integer to A/"(Ti,£), and <x% is the estimated 
uncertainty of the i-th TOA. The information matrix is defined as 

Lab = h -*r • (74) 

If only white noise is present in the TOA residuals, the inverse of the information matrix, M a f, = L~c, is the 
correlation matrix of the fitting parameters. The elements of the main diagonal of M a b give the variations 
of the measured parameters, (£„) = M aa , and the off-diagonal elements of M a b represent the correlations 



between them ( Bard 1974 ) 



Two kinds of BH-pulsar binaries seem likely to be present in our Galaxy. First, young binary pulsars 
located in the Galactic plane, where the pulsar is the result of the second supernova explosion. Numerical 



simulations indicate, that such a pulsar is expected to be in eccentric wide (P& J> 10 days) orbits (Lipunov 



et al. 1994 ). Secondly, a millisecond pulsar captured by a BH, most likely located in a globular cluster 



( Kulkarni, Hut, fc McMillan 1993 ) -PI The aim of our numerical estimates is to understand which parameters 
one can expect to measure in these two different types of BH-pulsar binaries within a reasonable time span 
of observations (<; 20 years). Here we assume that all measured TOAs, tj, are of a similar quality, i.e. all 
are the same (cr, = cttoa)- Furthermore, there is at least one timing observation every month after an 
initial set of 10 timing observations along one full orbit which give the Keplerian parameters required for 
subsequent exploration of the binary system. We further assume that the pulsar is in a highly eccentric orbit 
(e = 0.8) with a 10M Q BH companion. 

Figure 4 shows the estimated time span of regular timing observations, T b s , needed for the measurement 
of certain astrophysical parameters or effects, as a function of the orbital period, P^. As astrophysical 
parameters and effects we consider: 



9 The ratio of orbital binding energy to potential energy in the gravitational field of the cluster is only a very weak function 
of the companion mass, once m c > 1.4Mq. Thus, BH-pulsar binaries should have roughly the same probability to remain in or 
near the cluster after formation as NS-NS binaries, like PSR B2127+11C. 
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Total mass of the binary system — measurement of Cj. 
Mass of the BH — measurement of Co and 7. 

Testing the emission of gravitational waves — measurement of u, 7 and p, (cf. Taylor fc Weisberg 



1989| ). 



• Frame-dragging effect — measurement of x. 

• Spin of the BH (magnitude and orientation) — measurement of Co, 7, P,, x, and x. The mass deter- 
mination has to be done using 7 and p,. The accuracy for the total mass has to better than ~ eCos/Co, 
where s denotes the aspired accuracy for the BH spin measurement. 

In figure 4a we have assumed a timing accuracy of otoa ~ 100 fj,s which represents the optimum one can 
expect for timing observations of a young pulsar, with a typical rotational period of 100 ms to a few seconds. 
In figure 4b we have assumed the presence of a millisecond pulsar, i.e. otoa ~ 1 A<s. 

We conclude from figure 4a that for a young pulsar only in systems with orbital periods below 0.2 days 
the measurement of the spin of the BH (magnitude and orientation) is possible within a reasonable observing 
time span. For systems with orbital periods below 3 days one can expect to see the frame-dragging caused 
by the extreme Kerr BH, i.e. measure x. Systems with orbital periods below 10 days will allow a mass 
determination for the BH with an accuracy of better than 5% after less than 20 years of observing. The 
periastron advance will be a prominent effect even for orbital periods exceeding one year and thus the total 
mass of the binary system should be a measureable quantity for all orbital periods <1000 days. 

In the case of a millisecond pulsar in orbit with a 10 solar mass extreme Kerr BH, where one expects 
timing accuracies two orders of magnitude better than in the previous case, the corresponding time spans for 
observations are clearly shorter (see figure 4b) . Therefore the spin of the BH will be measurable for systems 
with orbital periods up to one day, and the frame-dragging will be seen even in systems with 10 days orbital 
period. For systems with orbital periods below a hundred days a good mass estimates for the BH companion 
will be possible, and the total mass of the binary system will be measurable even in very wide systems with 
orbital periods of a few years. 

The previous calculations were done for a BH mass of 10M Q , which we assume to be typical for stellar- 
mass BHs. In a 30A/q BH-pulsar binary spin-related effects are expected to be an order of magnitude 
larger, since the spin of a BH scales with its mass squared (see equation (H)). On the other hand, there 



might be a fairly large population of BHs in our Galaxy with masses slightly above l.5M@ (Brown & Bcthc 



1994). Our calculations show that only for very short orbital periods (Pf, < 0.2 days) and high orbital 
eccentricities, (e J> 0.8) one can expect to see the relativistic spin-orbit precession (is) within less than 
20 years of observing, given an uncertainty in the TOAs which is not worse than l/is. In this case the 
measurement of is, which gives a lower limit for Xi might have the potential to distinguish between a NS 



and a BH, since it seems that x ^ 0.7 for NSs of comparable mass (Laarakkers & Poisson 1997). It would 
require that the angle 9 is close to 90° (the spin of the BH lies almost in the orbital plane) and $0 is either 
close to 90° or close to 270°. Then, if \ is close to one for the BH companion, the minimum value for \ 
derived from the x measurement can exceed the maximum value for NSs. 

Let us point out that on time scales of several years, timing observations are usually contaminated 
by the increasing amount of red (low frequency) timing noise flCordes fc Helfand 1980j [Taylor 199l| , |Cy7e] 



1996, Kopcikin 1997b). If the red noise dominates in TOAs residuals, it makes the problem of measuring 



relativistic parameters more difficult. In particular for young pulsars in wide orbits, the separation of orbital 
effects and timing noise may cause a severe problem (cf. Wex et al. 1998| ). 



7. Measuring the quadrupole of the BH companion 

Equation (Q) implies that with the determination of the mass and the spin of a Kerr BH we also know 
its quadrupole moment. Therefore, if we are able to extract independently the quadrupole moment of the 
companion from our timing observations we could actually test whether the observed pulsar is orbiting a 
Kerr BH or another compact relativistic object like neutron or boson star. As discussed in Section 4 the 
quadrupole moment of the BH companion will lead to an additional precession in the angles $ and iff and, 
thus, to a secular time evolution in the parameters u> and x. Unfortunately these secular changes in the 
orientation of the orbit caused by the quadrupole moment of a BH companion are typically three orders 



of magnitude smaller than the changes caused by the relativistic spin-orbit coupling (see equation (49)). 
Thus, the chance of extracting the quadrupole moment of the BH from the measurement of the parameters 
to and x is unlikely. On the other hand, the anisotropic nature of the quadrupole component of the external 
gravitational field (||) will lead to characteristic short-term periodic effects every time the pulsar gets close to 
the oblate companion. A detailed investigation of these short-term effects can be found in Wex (1998). These 
short-term periodic effects lead to a unique signature, which will show up in the post-fit TOA residuals if one 
uses a timing model that accounts only for the secular changes of the precession and the short-term periodic 



effects caused by the mass-monopole (Damour k, Deruelle 1986) and spin (Wex 1995) of the BH. Figure 5a 



shows the result of numerical simulations of TOA residuals for a pulsar orbiting a 30 M Q BH in an 0.1 day 
orbit with an eccentricity of 0.9. This represents certainly an unrealistic case from the observational point 
of view since the lifetime of such a system due to gravitational radiation damping is just 10 4 years which 
makes a discovery of such a system extremely unlikely. In addition, the accuracy needed to measure this 
effect would require a nano-second timing precision. Even for the best millisecond pulsars the present timing 
precision is still more than two orders of magnitude worse and, thus, we conclude that the quadrupole of a 
stellar-mass BH companion is not measurable in timing observations. However, since the quadrupole moment 
of a BH scales with the mass cubed (see equation ^[), these short-term periodic effects under discussion gain 
importance for very massive BHs which could sit in the center of globular clusters as a result of a collapse 
of the cluster core. Figure 5b gives numerical simulations of TOA residuals for a pulsar orbiting a 10 4 M & 
BH in a 10 day orbit with an eccentricity of 0.9. A timing accuracy of a few hundred nanoseconds would 
allow the measurement of the quadrupole moment. Of course, even higher BH masses and/or smaller pulsar 
orbits will make the quadrupole a prominent feature in the TOAs. 



8. Summary 

In this paper we have examined how timing observations of a pulsar orbiting a rotating stellar mass BH 
can be used to study the physical properties of the BH companion. We pointed out that the measurement of 
two post-Keplerian parameters can lead to a mass determination of the black hole with very high accuracy, 
far better than it is possible for present black hole candidates. It was shown that the frame dragging 
propagation effect discussed by Laguna & Wolszczan (1997) will not be separable from the bending delay for 
those BH companions less than a hundred solar masses, since it would require determining the orientation of 
the pulsar's spin with an accuracy that is not likely possible with pulse structure analysis. We gave detailed 
calculations about the (secular) influence of the BH rotation on the orbital motion of the pulsar. In particular 
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the long term behavior of this relativistic spin-orbit coupling can be described by two angles, defined with 
respect to the invariable plane, which change linearly in time. For an observer this linear-in-time precession 
converts into a non-linear in time evolution of the orbital inclination and the longitude of periastron. In a 
timing model for such a binary system, this can be taken into account by fitting for the parameters x, x, u>, 
and u> as is done for pulsars with main-sequence star companions where the classical spin-orbit coupling is 



the cause of the precession (see Wcx 199S). The measurement of x, x, and d) could, in principle, lead to the 
determination of the BH spin, if the masses of the BH and the pulsar can be determined without making use 
of the advance of periastron, using, for example, the measurement of the Einstein delay (7) and the orbital 
period decay (Pb)- In order to do so, a fractional measurement precision of <J 10 -4 is required for u>, 7, and 
Pb- If this measurement precision is not achieved, or if there are Galactic and kinematic effects present which 
are not sufficiently understood, then only the projection of the BH spin onto the orbital angular momentum 
will be measurable. If the timing data allow just the determination of x then we can calculate only a lower 
limit for the projected BH spin, which, of course, is also a lower limit for the BH spin. However if we have 
just the lower limit for the spin, the presence of the x would already indicate frame dragging if an oblate 
non-compact companion can be ruled out (for example by optical observations, the absence of eclipses of the 
pulsar, absence of tidal effects, etc.). On the other hand, the method under discussion can be used to rule 
out a Kerr BH companion as has been demonstrated here in the case of PSR J0045-7319. 

We have studied the measurability of the relativistic spin orbit coupling for two different cases. First, a 
young pulsar orbiting a IOMq extreme Kerr BH and second, a millisecond pulsar orbiting a IOMq extreme 
Kerr BH. In the first case only for binary pulsars with orbital periods less than three days one can expect 
to measure x, i.e. the frame dragging caused by the rotating BH, with sufficient accuracy, but in the second 
case a measurement of x seems very likely. For BH-millisecond pulsar binaries with an orbital period below 
one day one can expect a measurement of x (and u>) within a reasonable time span of observations, and one 
will be able to fully determine the BH spin parameters (magnitude and orientation). However, if the pulsar 
is located close to the core of a globular cluster one might not have the necessary understanding for the 
gravitational cluster potential in order to take properly into account its contribution to the measured binary 
parameters. For a 1.5M BH-millisecond pulsar binary with Pb <^ 0.2 days and e J> 0.8 one might be able 
to measure is- In principle, this could help to exclude a neutron star, if the magnitude and orientation of 
the BH companion is such that \ sin 9 sin $0 i> 0.7, since NSs seem to have \ <; 0.7. 

Finally we pointed out that the quadrupole field of the rotating BH will lead to a distinctive signature 
in the post-fit TOA residuals during each periastron passage. However, our numerical simulations suggest 
that this would be observable only for very massive BHs, typically more than 10 4 M Q . Whether BHs of this 
size do exist, e.g. in the center of globular clusters, is still unclear. The discovery of a pulsar orbiting the 
super-massive BH in the center of our Galaxy could provide a possible setting for these measurements. 
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Fig. 1. — Mass ranges of companion stars in pulsar binaries as estimated from timing observations. The 
solid circles give the most likely value for the companion mass while the error bars indicate a 95% confidence 
level. The lower horizontal dashed line indicates the Chandrasekhar mass limit, the upper one indicates 
the 3 Mq upper mass limit for neutron stars. There are only five binary pulsars where the plotted mass 
range of the companion exceeds 3 Mq. The first two binary pulsars from the left (B0655+64, J1022+1001) 
have optically detected white dwarf companions (Kulkarni 1986; Lundgren, Foster, & Camilo 1996). The 
third binary pulsar (J0045— 7319) has a 9 Mq B-star companion (Bell et al. 1995), and the fifth binary 
pulsar (B1259-63) is in orbit with a 10 Mq Be-star (Johnston et al. 1994). Only the nature of number four 
(B1820— 11) is unclear. Mass functions and masses were taken from (Deich & Kulkarni 1996, Nice, Sayer & 
Taylor 1996, Stairs et al. 1998, Kaspi, Taylor & Ryba 1994, Taylor et al. 1993, 1995). We further assumed 
that the mass of the pulsar is above 1.3 Mq 
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Fig. 2. — Orientation of the BH-PSR system in the sky. The vector S is the spin of the BH, R is the vector 
pointing from the BH to the pulsar. The angle A. is that between the BH spin and the direction of the 
line-of-sight given by the unit vector K . The angle r], is that between the projection of the BH spin on to 
the plane of the sky and the ascending node of the binary orbit. 



plane of the ski) 



Fig. 3. — Definition of angles in the total- angular- momentum reference frame. The invariable (X-Y) plane 
is perpendicular to the total angular momentum J = L + S. The line-of-sight vector Ko is in the Y-Z plane. 
The vector J is a conserved quantity and, if averaged over one full orbital period, the absolute values |L| and 
S are constant. Thus, the angles ij, Oj, and 9 are fixed. The angles <J> and ^ change linearly with time. 
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Fig. 4. — Time span of observations needed (after 10 timing observations along a full orbit) to determine the 
total mass of the BH-PSR system with 5% accuracy (full), the mass of the BH with 5% accuracy (dashed), 
the emission of gravitational waves with 1% accuracy (dot-dash), the precession caused by frame dragging 
with 5% accuracy (dotted), and the spin of the BH with 5% accuracy (dash-dot-dot-dot) as a functon of the 
orbital period P&. Estimations were done for a 10 solar mass extreme Kerr BH with 9 = 45° and $o = 45°, 
an orbital eccentricity of 0.8 and a timing accuracy, otoa, 01 100 ^s (upper figure) and 1 [is (lower figure). 
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Fig. 5. — Upper figure: Typical signature in the timing residuals caused by the quadrupole moment of a 30 
Mq BH companion. We used Pb = 0.1 days and e = 0.9. The inclination of the BH spin with respect to the 
orbital plane (the angle 0) was assumed to be 70°. Lower figure: Timing residuals caused by the quadrupole 
moment of a 10 4 Mq BH companion. We used Pb = 10 days and e = 0.9. The inclination of the BH spin 
with respect to the orbital plane (the angle 9) was assumed to be 70°. 



